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Abstract 

The top quark forward-backward asymmetry A‘,,, measured at the Tevatron is above the Stan- 
dard Model prediction by more than Ze deviation, which might be a harbinger for new physics. 
In this work we examine the contribution to A‘,, in two different new physics models: one is 
the minimal supersymmetric model without R-parity (RPV-MSSM) which contributes to A^. via 
sparticle-mediated t-channel process dd — tt; the other is the third-generation enhanced left-right 
model (LR model) which contributes to Atp via Z'-mediated t-channel or s-channel processes. 
We find that in the parameter space allowed by the tt production rate and the tt invariant mass 
distribution at the Tevatron, the LR model can enhance A‘,, to within the 2c region of the Teva- 
tron data for the major part of the parameter space, and in optimal case A‘,, can reach 12% 
which is slightly below the lo lower bound. For the RPV-MSSM, only in a narrow part of the 
parameter space can the A" couplings enhance Ar to within the 20 region while the Al couplings 


just produce negative contributions to worsen the fit. 


PACS numbers: 14.65.Ha,14.80.Ly,11.30.Hv 


I. INTRODUCTION 


As the heaviest fermion with a mass at weak scale, top quark is speculated to be a sensitive 
probe of new physics beyond the Standard Model (SM) il The precision measurement of its 
properties, now being performed at the Tevatron and will be continued at the CERN Large 
Hadron Collider (LHC), will either unravel or further constrain the new physics related to 
the top quark. 

So far the production rates of top pair and single top measured at the Tevatron are in 


good agreement with the SM predictions, but a more than Ze deviation is reported in the 


forward-backward asymmetry A‘, in top pair production, which is defined by 
FB 


APP. — N,(cos@ > 0) — N;(cos0 < 0) (1) 
FB N,(cos@ > 0) + N;(cos0 < 0)’ 


with 0 being the angle between the reconstructed top quark momentum and the proton beam 
direction in tt rest frame. In the SM, this asymmetry gets dominant contribution from the 
next-to-leading-order QCD correction and was found to be several percent: Afp(SM) = 
5.0 41.5% [2]. Compared with its experimental value Atp (exp) = 19.3 + 6.996 measured by 
the CDF and D0 collaborations B. this SM prediction is below the measured value by more 
than 2c deviation. 

Such a discrepancy might be a new physics footprint in top quark sector and has been 
studied in several new physics models, where the Kaluza-Klein excitations in extra dimen- 
sions HI. the presence of new gauge bosons (Z', W’, axigluon) HE or new scalars H are 
utilized to try to explain the discrepancy. Noting that the mechanisms proposed in these 
literatures can also be realized in some popular new physics models, we in this work study 
the asymmetry A‘,, in the supersymmetric models and the left-right models 

For the minimal supersymmetric standard model (MSSM) [9], the SUSY influence on tt 
production comes from loop effects, in which the SUSY-QCD effects are dominant over the 
SUSY-EW effects . But among the SUSY-QCD one-loop diagrams only the box diagrams 
contribute to the asymmetry A‘, and, consequently, the contribution is negligibly small (see 
the discussion in |4]). Therefore, in our analysis we consider the R-parity violating MSSM 
(RPV-MSSM) which allows for tree-level contribution from t-channel process dd — tt by 
exchanging a color-singlet slepton or a color-triplet squark. For the general left-right models, 


since the predicted new gauge bosons are usually at TeV scale and unlikely to affect A‘, 


significantly, we here consider a special left-right model called the third-generation enhanced 
left-right model l3. This model predicts a new gauge boson Z’ which contributes to Ain 
via t-channel or s-channel processes. 

'This paper is organized as follows. In Sec. II we describe the calculation of the asymmetry 
Abp in the RPV-MSSM and present some numerical results and discussions. In Sec. III 
we perform similar analysis in the third-generation enhanced left-right model. Finally, some 


discussions and the conclusion are presented in Sec. IV. 


II. Aj IN R-PARITY VIOLATING MSSM 


In the popular MSSM, the invariance of R-parity, defined by R = (—1)?5*?P*7 for a field 
with spin S, baryon-number B and lepton-number L, is often imposed on the Lagrangian 
in order to maintain the separate conservation of B and L. Although R-parity plays a 
beautiful role in the phenomenology of the MSSM (e.g., forbid proton decay and ensure a 
perfect candidate for cosmic dark matter), it is, however, not dictated by any fundamental 
principle such as gauge invariance and there is no compelling theoretical motivation for it. 
The most general superpotential of the MSSM consistent with the SM gauge symmetry and 
supersymmetry contains R-violating interactions which are given by L 


1 1 
Wg = ^ui; Ex + AjgyLiQ;Dy + hin" H Djs Di. + Bio Hs, (2) 


where 4, j,k are generation indices, c denotes charge conjugation, a, 8 and y are the color 
indices with e?" being the total antisymmetric tensor, Hə is the Higgs-doublet chiral super- 
field, and L;(Q;) and E;(U;, D;) are the left-handed lepton (quark) doublet and right-handed 


lepton (quark) singlet chiral superfields. The dimensionless coefficients A;;; (antisymmetric 


in ¿ and j) and X;,, in the superpotential are L-violating couplings, while Ai, (antisymmetric 
in j and k) are B-violating couplings. 

The expression of Wy implies that both A7;,and Xj, can induce new top quark interac- 
tions. In terms of the four-component Dirac notation, the interactions involved in top pair 


production are 
"ü gr. 3 1 7i age Fj*—i c 
L= Aldin, E adr uud, + dudi] + h.c., (3) 


and these interactions can contribute to the forward-backward asymmetry by the diagrams 


d(p2) (pa) d(p) t(pa) 
FIG. 1: Feynman diagrams contributing to A in the RPV-MSSM, where i and dk, denotes a 


left-handed slepton in i-th generation and right-handed squark in k-th generation, respectively. 


TABLE I: The upper bounds on the couplings Aa (i = 1,2,3) and Aan, (k = 2,3) l4. 


couplings bounds sources 
Aen 0.03 m; ; /(100 GeV) Ou: (C's) 
Agi 0.18 mj. /(100 GeV) Vig 
Aen 0.26 mg. /(100 GeV) K > nvv 
Asik 0.97 mg, /(100 GeV) RZ 
Asik 1.25 perturbativity 


shown in Fig[] The corresponding amplitudes are then given by 
pau ) Pau(d)v(d) Prv(t) (4) 
(pı — p3)? — mis l 
(t) Yu Prv(t)o(d)y" Pru(d) 


: u 
= eeh eo hann) 2l(p — pa)? — m? | ? p 
dkR 


MA yaly = el äel Aen 


MEP 


dd—1tt VIV 


with a, 8, p, c and A being color indices of the quarks and squarks. These amplitudes affects 
Atp by interfering with the QCD amplitude dd — g* — tt and also by its own square. In 
our results presented below, we have included the SM contribution to Ae and considered 
only one coupling non-zero each time. 

The SUSY parameters involved in the calculation are the couplings M4, and Ana, as well 
as sparticle masses. So far both theorists and experimentalists have intensively studied the 
phenomenology of these couplings in various processes and obtained some bounds |14]. 
In Table [I] we list the relevant bounds and, as can be seen, in case of heavy squarks, these 
bounds are quite weak. Note that for \3,, we do not use the stringent bound from n — ñ 
oscillation because they depend on additional SUSY parameters which are not involved 
in our processes. 


In Fig2] we show the dependence of Ac on the relevant sparticle masses. Here we 
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FIG. 2: Ar versus the sparticle mass in the R-parity violating MSSM. The left (right) panel is 


induced by ai pu ,) which are fixed at their maximally allowed values. The SM contribution to 


Ae is also included. 


assume all squarks and sleptons are degenerate and sum over the contributions from different 
generations. The couplings Aa and \3,, are fixed at their maximally allowed values which, 


as shown in Table I, vary with squark mass. The SM parameters are taken as [15 
m, = 172.5 GeV, mz = 91.19 GeV, sin Ow = 0.2228, o,(m;) = 0.1095, a = 1/128. (6) 


Fig[2]indicates that both A4; and Au can give a negative contribution to Afp and in 
the worst cases, they decrease Air by 2% and 5% respectively, which are comparable in size 
with 44.2", In a narrow part of the squark mass (250 ~ 400 GeV), the coupling A3,, can 
also give a positive contribution to enhance A55 to within the Ze region of the Tevatron 
data. 

We note that the scalar-mediated contributions have been analyzed in a model indepen- 
dent way in [7| and the results were presented for color-singlet, -triplet, -sextet and -octet 
scalar respectively. We checked that our analytic results are in agreement with for the 
color-singlet and -triplet cases except that our study is restricted in a specified model, the 


RPV-MSSM. 


Before ending this section, we give a comment on A‘,p in the top-color assisted technicolor 
model (TC2) Si This model predicts some composite bosons, 7? and z 7. with mass at weak 


scale and having large Yukawa couplings to top quark. As a result, Ar gets additional 
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contribution from t-channel processes uu(cc) — tt by exchanging a color-singlet scalar 79, 
which is similar to Fig[1](a) in RPV-MSSM. Noting the up quark content in proton is larger 
than the down quark content, one may expect a larger effect on Air in TC2 model than 
in RPV-MSSM. This is not true because for uŭ — tt in TC2 model, its contribution is 
proportional to ug — tr mixing which is determined by the triangular texture of the up-type 
quark mass matrix and turns out to be very small . As to cc > tt, although cg —1g mixing 
may be sizable A it actually gives no contribution to An because the distributions of c 
and c in proton are approximately same and so the initial state of this process is symmetric 
under the exchange of c and c. In fact, we numerically calculated all TC2 contributions 
to tt production at the Tevatron, which include the s-channel processes gg, bb > zb — tt, 
t-channel processes bb — tt by exchanging 7; and uü,cc — tt by exchanging 7°, and we 


found that TC2 model can change oj by at most 200 fb and A‘,, by order of 10 *. 


II. Abp IN THE THIRD-GENERATION ENHANCED LEFT-RIGHT MODEL 


In the third-generation enhanced left-right model 13. the gauge group is SU(3)o x 
SU(2), x SU(2)n x U(1)pg.r with gauge couplings gs, gr, gr and g, respectively. This 
model differs from other left-right models by having right-handed gauge bosons couple pre- 
dominantly to the third generation fermions. Noting that the right-handed gauge boson Zp 
will mix with the standard Zp to form mass eigenstates Z and Z', one can write down the 
neutral gauge interactions of quarks as 


gL ' 
ll (9v — gays )qlcos £z Zp — sin £z Z) 
Ow 


Ly = — 
e 2 cos 


1 4 2-2 . 
+ > tan On(saLY"aL + gni" uni — 34 ni" dni)(sin £z Z, + cos £z Z,,) 


— + (tan dr + cot Or) (üni^" Vg Vg;uRj — dria” Vis; Vg; dry) (sin £z Z, + cos zZ) (7) 
where tan fr = g/gr, gy = 9 608 OR = grsin Op, £z is the mixing angle between Zp and Zo, 
and VM are unitary matrices which rotate the right-handed quarks ug; and dg; from inter- 
action basis to mass eigenstates. This model also predicts new charged gauge interactions 
of quarks, but since their effect on A‘,, is small, we do not consider them here. Note that 
in Eq.(1]), q and qz are summed over all quarks, and the repeated generation indices 7 and 
j are also summed. 


Eq. indicates that the Z'u;u; interaction is large when gg >> gy or cot@p > 1. 
j 
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u(pi) t(ps) u(pi) t(ps) u(pi) t(p3) u(pi) t(p3) 
Z Z 
U(p2) (a) t(pa) lpo) (b) t(pa) ü(pa) (o) t(p,) ü(p3) (d) t(pa) 


FIG. 3: Feynman diagrams contributing to 44%, in the third-generation enhanced left-right model. 
FB 


This feature can be utilized to enhance An by the t-channel process uŭ — tt if (Vg), is 
moderately large. In PI V£ and V7 are assumed to be nearly diagonal to satisfy the severe 
flavor-changing neutral-current (FCNC) constraints. However, as pointed out in H a sizable 
Uug-—tg mixing does not conflict with these constraints given that the other flavor mixings are 
suppressed. Here we further point out that this pattern of flavor mixings does not necessarily 
require the up-top element in up-type quark mass matrix M,, to be much larger than other 
off-diagonal elements. For example, assuming (Vz),; = 0.2, (Vg)a = 0 and (V#)ue = 0, we 


with M? 


numerically solve the equation VŠ MIM,VE = M? diag 


— 2 NU 2 
diag = diag(m£,mi,m;, and 


we find it possible that (M), is several times larger than (ML, Ju. 
In the third-generation enhanced left-right model, beside the dominant QCD contribution 
qq — g* — tt, diagrams shown in Figl3] also contribute to tt production at the Tevatron, 


and the amplitudes of these new contributions are 


E e V'ü(t)ouloz, Pr + gg PRjv(t)v(u)y" [gg Pu + ggg Pn]u(u) 
Ma = Sat uns, Pe eher D rëselt ` e 
P (p1 + p2)? — mz 
My A.T E chen + dy pPa OR oy Pi + 9% T dO NN 
GEES (pi + p2)? m?, iD2mg i 
2 _ E 
€ u vu 2U(t)y, Pgu(u)v(u)y" Pnv(t 
Me = Wap9g0 (=) [£z s, (cot Op + tan 0g)]^ [Vil Va ee (io) 
Ww =. m Z 
2 ` ` 
e 2 wen yu Ut) Yu PRu(u)v(u)Y" Pgv(t) 
Ma = ióopógs | ——— w 0 ———————É——-—, (M 
"et ( 5 De (cot e+ tan ba) Vi viu P DD R, (am 


where Uz is the Z’ width obtained by adding all decay channels of Z', Sw = sin Hu, Cy = 


. : tu t, t t, 
cos Ôw, and the coupling coefficients 977, däm, 971, and 92 are defined as 


1 
Iz, —1- zou — g$v tan g£z, (12) 
4 4 
IZR = —39w T gv tan frt (13) 
4 1 
dou = E = gs, tan OREZ + Sw cot OrREz, (14) 


4 1 
97, = (1— ssn )éz + zë tan Op, (15) 


2 3 
4 4 
gg ge? -+ 39v tan Op, (16) 
4 1 
= ës? + gv tan Üp£z — Sy cot Op. (17) 


Note that among the four amplitudes, only M. and Mj interfere with the QCD amplitude. 

In this model the new parameters €,, cot 0g, Mz and (Vg), are involved in our cal- 
culation. Constraints on these parameters were discussed in i and it was found that 
0 < & € 0.02 and cot0g < 20. As for Mz, we should note that the constraints from 
CDF search for Z' and from the global fitting of the electroweak precision data are 
invalid here since these constraints arise mostly from the processes involving the first- or 
second-generation of fermions. So far the pertinent bound comes from ee" — bb at LEP-II, 
which requires Mz = 460 GeV for cot fr > 10 [12]. In following analysis, without stating 


explicitly, we include the QCD contribution to Air and scan the new free parameters in the 


following ranges 


500 GeV € Mz < 2000 GeV, 0< z €0.02, 10x cot0g € 20, 0.1 € (Vg). € 0.2 
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FIG. 4: Scatter plots of A‘, versus cot Ôp and m — for the no-mixing case in the LR model. 


In our discussion we consider two cases, i.e., with and without ug — tr mixing. For the 
no-mixing case, new contributions to An only come from diagrams (a) and (b) of Figl[3l 


and the dependence of A‘,, on cot Og and mz are shown in Figli] This figure shows that a 
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light Z’ with relatively small cot fo can enhance Aj,, to the 26 region of its experimental 
value. This is because small cot Og can enhance Z'üu interactions, and a light Z’ can make 
the resonance effect of the diagram Fig[3[b) on tt production more significant. FigH also 
shows that in the no-mixing case, the Z’ contribution to Ae can only reach 1%. This is 
due to the smallness of Z'uu couplings (see Eq.(14) and Eq.(15)) and the absence of the 
interference of diagram (a) and (b) in Fig] with the dominant QCD amplitude. 


Without the Constraints With the Constraints 
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FB 


A 


12 14 16 18 10 | 12 14 16 | 18 20 
coto, cote, 


FIG. 5: Scatter plots of An versus cot 0g in the mixing case of the LR model. The left (right) 


panel is without (with) considering the constraints from our and Mz. 


Without the Constraints With the Constraints 
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FIG. 6: Same as Fig. [5] but projected in Ac versus Mz, plane. 


In Figs[5land [6] we show Air versus cot 0g and Mz, respectively in the mixing case. The 
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left panels of these figures show that the effect of Fig[3] (mainly diagram (d)) is potentially 
very large, pushing the value of Atp up to 50%. In this case, the effects on the total tt 
production rate and the tt invariant mass distribution are also large and so it is necessary 
to consider constraints from these observables measured at the Tevatron. Now based on the 
CDF 4.6 fb^! luminosity data, the measured total cross section is a? = 7.50 + 0.31,44 + 
0.34,,,; = 0.154, pb for m, = 172.5 GeV bi Combining errors in quadrature, one can get 
cz" = 7.50 + 0.48 pb, which is in good agreement with the SM prediction ou = 7.519? 
pb [21]. The invariant mass distribution was also measured by CDF, and the results are 
presented in nine bins of Mẹ |22|. When we calculate these observables, we multiply an 
overall K-factor 1.329 for the dominant tree-level QCD contribution bi. The right panels 
of Figs[5] and [6] are then obtained by requiring the total cross section and the differential 
cross section in each bin to be within the 2e regions of their experimental values. These 
results shows that even with the constraints, the third-generation enhanced left-right model 
can still enhance Aj, to 12% (well above the 2ø lower bound), but it can not enhance Abp 
to within the Lo region of its experimental value. Note that unlike the no-mixing case, large 
contribution to An comes from the region where cot Oz is large. The reason is that in the 
mixing case, dominant contribution arises from diagram (d) of Fig] and this contribution 


is proportional to (cot 05 + tan OR)’. 


IV. DISCUSSION AND CONCLUSION 


From our study and the previous works H nn A we can learn that to make sizable 
contribution to Atp, the following two conditions must be satisfied. One is the initial state 
must be uŭ and/or dd and in case that only dd initiated contribution is responsible to explain 
A*. 5, the involved interaction must be strong enough to compensate the suppression of the 
parton distribution of down quark in proton. The other is the amplitude of the tt production 
must contain terms proportional to cos 0 with 0 being the angle between the reconstructed 
top quark momentum and the proton beam direction in tt rest frame, or in other words, 
contain terms like (p, - pi) (pa - pi) — (pa: Pt)(Pu - pi). This requirement implies that new 
physics affect tt production through the following ways: 


(1) Through s-channel process qq — tt by exchanging a gauge boson H H If this process 


interferes with the the s-channel QCD process qd — g* — tt, the interaction of 
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the gauge boson with q and t must have axi-vectorial component. Examples in this 
direction are the presence of Kaluza-Klein excitation of gluon in extra dimension |4] or 
the axigluon |5]. If this process does not interfere with the s-channel process, to affect 
Abp by itself, the interaction must have both vectorial and axi-vectorial component, 
like what we studied in diagram (b) of Fig[3] From yet known studies we can infer 
that it seems difficult for the latter case to enhance Atp significantly without spoiling 


the constraints from o, and Ma. 


(2) Through t-channel process gg — tt by exchanging a vector boson or a scalar ld H 
which interferes with the QCD process qq — g* — tt. In this way, scalar is less 
efficiency than vector boson in explaining Air given that they have the same coupling 
strength and mass. This is because for the scalar case, there is a competition between 
spin correlation and the Rutherfold singularity Të. Moreover, as shown in [7], the 
scalar-mediated contributions can be categorized by the transformation property of 
the scalar under SU(3).. If the scalar is color-triplet or -sextet, there exists large 
parameter region to explain Atp within 1o and at the same time to remain oe within 
the experimental errors, while for color-singlet or -octet scalar, it is very difficult to 
"c. a large positive contribution to Ak without spoiling the constraint from ou 


. In our work, we checked this conclusion for color-singlet and -triplet cases. 


From our study, we can also learn that, although the top quark pair productions at the 
Tevatron and the LHC may be sensitive to new physics, the effects of new physics (like the 
popular MSSM or TC2 models) are usually not so large to be well above the experimental and 
theoretical uncertainties . À complementary or even more sensitive probe for new physics 
effects in top quark sector is top quark FCNC processes, which are extremely suppressed and 
unaccessible in the SM but can be greatly enhanced by several orders to reach the observable 
level in some new physics models like the MSSM SÉ the TC2 model ba. 

In summary, we in this work calculated the new physics contribution to the top quark 
forward-backward asymmetry Ar at the Tevatron in two different models: the minimal 
supersymmetric model without R-parity (RPV-MSSM) and the third-generation enhanced 
left-right model (LR model). We found that in the parameter space allowed by the tt produc- 
tion rate and the tt invariant mass distribution at the Tevatron, the LR model can enhance 


A. to within the 20 region of the Tevatron data for the major part of the parameter space, 
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and in optimal case A can reach 12% which is slightly below the Le lower bound. For the 
RPV-MSSM, only in a narrow part of the parameter space can the A" couplings enhance 
A to within the 2c region while the A couplings just produce negative contributions to 
worsen the fit. Noting that the R-parity conserving interactions in the MSSM are unlikely to 
give large enough contribution to Ak, we conclude that the MSSM with (without) R-parity 
will be disfavored (favored in case of A" couplings) if the discrepancy of A‘,, persists with 
more forthcoming data accumulated by the Tevatron. We also checked the top-color assisted 
technicolor model and found that it gives negligibly small contributions to Abp and thus 


unlikely to explain the Tevatron data. 
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